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Abstract 

We investigate the transition formfactors for the B — ► K*l + l~(l = [i, r) decay in 
the light cone QCD. It is found that the light cone and 3-point QCD sum rules analyses 
for some of the formfactors for the decay B — * K*l + l~ lead to absolutely different 
q 2 dependence. The invariant dilepton mass distributions for the B — > K*fi + fi~ and 
B — ► K*t + t~ decays and final lepton longitudinal polarization asymmetry, which 
includes both short and long-distance contributions, are also calculated. 
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1 Introduction 



Experimental observation |J of the inclusive and exclusive radiative decays B — > X s 7 and 
B — > K*^ stimulated the study of rare B decays on a new footing. These flavor changing 
neutral current (FCNC) b — > s transitions in the SM do not occur at tree level and appear 
only at loop level. Therefore the study of these rare S-meson decays can provide a means of 
testing the detailed structure of the SM at the loop level. These decays are also very useful 
for extracting the values of the Cabibbo-Kobayashi-Maskava (CKM) matrix elements ||, 
as well as for establishing new physics beyond the Standard Model ||. 

Currently, the main interest on rare S-meson decays is focused on decays for which the 
SM predicts large branching ratios and can be potentially measurable in the near future. 
The rare B — > Kl + l~ and B — > K*l + l~ decays are such decays. The experimental situation 
for these decays is very promising [|J, with e + e~ and hadron colliders focusing only on 
the observation of exclusive modes with I = e, (/, and r final states, respectively. At 
quark level the process b — > sl + l~ takes place via electromagnetic and Z penguins and W 
box diagrams and are described by three independent Wilson coefficients C 7 , C 9 and C w . 
Investigations allow us to study different structures, described by the above mentioned 
Wilson coefficients. In the SM, the measurement of the forward-backward asymmetry 
and invariant dilepton mass distribution in b — > ql + l~ (q = s, d) provide information on 
the short distance contributions dominated by the top quark loops and are essential in 
separating the short distance FCNC process from the contributing long distance effects 
and also are very sensitive to the contributions from new physics ||. Recently it has been 
emphasized by Hewett J7| that the longitudinal lepton polarization, which is another parity 
violating observable, is also an important asymmetry and that the lepton polarization in 
b — > sl + l~ will be measurable with the high statistics available at the B-factories currently 
under construction. However, in calculating the Branching ratios and other observables in 
hadron level, i.e. for B — > K*l + l~ decay, we have the problem of computing the matrix 
element of the effective Hamiltonian, Tt e ff, between the states B and K* . But this problem 
is related to the non-perturbative sector of QCD. 

These matrix elements, in the framework of different approaches such as chiral theory 
||, three point QCD sum rules method ||, relativistic quark model by the light-front 
formalism |lCj, [XT], have been investigated. The aim of this work is the calculation of these 
matrix elements in light cone QCD sum rules method and to study the lepton polarization 
asymmetry for the exclusive B — > K*l + l~ decays. 

The effective Hamiltonian for the b — ► sl + l~ decay, including QCD corrections [12-14] 
can be written as 



AC 10 

which is evolved from the electroweak scale down to \x ~ mj, by the renormalization group 
equations. Here represent the relevant CKM matrix elements, and Oj are a complete 
set of renormalized dimension 5 and 6 operators involving light fields which govern the 
b — > s transitions and Cj(/i) are the Wilson coefficients for the corresponding operators. 
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The explicit forms of Cj(/i) and Oi(/x) can be found in [12-14]. For b 
effective Hamiltonian leads to the matrix element 



sl + l decay, this 



M 



GfOL 



7T 



v tb v; s 



c e 9 ff s Ll ^b L m + c w s L ^b L h^ 5 i 



C 7 _ - 
2—sia IMU q u (m b R + m s L)b l^l 
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where q 2 is the invariant dilepton mass, and L(R) = [1 — (+)7s] /2 are the projection 
operators. The coefficient Cg (//, q 2 ) = Cg(fi) + Y(fi, q 2 ), where the function Y contains 
the contributions from the one loop matrix element of the four-quark operators and can be 
found in [12-14]. Note that the function Y(fi, q 2 ) contains both real and imaginary parts 
(the imaginary part arises when the c-quark in the loop is on the mass shell). 

The B — > K*l + l~ decay also receives large long distance contributions from the cascade 
process B — > K*ij)V> — > K*l + l~ . These contributions are taken into account by introducing 
a Breit-Wigner form of the resonance propogator and this procedure leads to an additional 
contribution to Cg of the form [I15J 



2tt 



a' 



E 



V=ip, ip' 



m v T(V -> l+r) 
, (q 2 - m\r) - im v T x 



As we noted earlier, for the calculation of the branching ratios for the exclusive B — > 
K*l + l~ decays, first of all, we must calculate the matrix elements (i^*|s7 M (l — 7s)g|-B) 
and (K* \sia ^p" (1 + / -f 5 )q\B). These matrix elements can be parametrized in terms of the 
formfactors as follows (see also 0): 
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where e* is the polarization vector of K*, p + q and p are the momentum of B and X* 
and P^ = (2p + q) . The formfactor A3 (g 2 ) can be written as a linear combination of the 
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formfactors A\ and A 2 (see 



2 m B + m K » 2 m B - m K * 2 
Mq } = 2m K . Mq } 2^^ Mq } ' (5) 

with the condition A^(0) = Aq(0). In calculating these formfactors we employ the light 
cone QCD sum rules. 



2 QCD Sum Rules for Formfactors 

According to the QCD sum rules ideology, in order to calculate the formfactors we start by 
considering the representation of a suitable correlator function in terms of hadron language 
and quark-gluon language. Equating these representations we get the sum rules. For this 
purpose we choose the following correlators. 



1 / d 4 xe i,lx {K*(p)\s(x) % (l - l5 )b(x)b(0)i l5 q(0)\0) , (6) 



1 / d 4 xe t(lx (K*(p)\s(x)ia^(l + 75 )6(x)&(0)z 75 g(0)|0) . (7) 



ng°(p,g) 

Here the first correlator is relevant for the calculation of the formfactors V(q 2 ), Ai(q 2 ), A 2 {q 2 
and A (q 2 ) and the second one for Ti, T 2 and T 3 . 

The main task in QCD is the calculation of the correlation functions (6) and (7). This 
problem can be solved in the deep Euclidean region, where both virtualities q 2 and (p + q) 2 
are large and negative. The virtuality of the heavy quark in the correlators (6) and (7) is 
large, of order m| — (p + q) 2 , and one can use the perturbative expansion of its propagator in 
the external field of slowly varying fluctuations inside the vector meson. Then, the leading 
contribution is 



n«(p,g) 



d 4 xd 4 k e i( - q - k)x 



(2tt) 4 (m 2 -k 2 ) 
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(8) 



{2K)i f/;/2 // 2 ) 9 l '(^|s(x)^(l + 75)(^ + m 6 )75g(0)|0> . (9) 

It is obvious from the above expressions that the problem is reduced to the calculation of 
the matrix elements of the gauge-invariant non-local operators, sandwiched in between the 
vacuum and the meson states. These matrix elements define the vector meson light cone 
wave functions. Following [TR [T7[ we define the meson wave functions as: 



(10) 
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The functions <fi±(u, /i 2 ) and 4>\\(u, fi 2 ) give the leading twist distributions in the fraction 
of total momentum carried by the quark in the transversaly and longitudinally polarized 
meson, respectively In (T7j it was shown that 



9l{u) 



l + (2w-l) 2 J , (13) 
6u(l-tt), (14) 

which we use in the numerical analysis. For the explicit form of <f)±(u, fi 2 ) we shall use the 



9l{u) 



results of [117 
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(16) 



Here b = yiVc — jR/, and 



n+1 i 

7» = ^|1+4E7 

3=2 J 



where C F = 



As in [|T7j, we will use the following values for the parameters appearing in eqs.(10)-(12) 
and eq.(15) : 



= 210 MeV, af*(/i = m b ) = 0.57, of *(// = m b ) = -1.35 
and tig (/i = m&) = 0.46 , 
4>\\(u, fi 2 ) = 6u(l — u). 
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Using eqs. (10-12), we get the following results from eq.(8) and eq.(9) for the theoretical 
part of the sum rules: 
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where 



and 
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Let us turn our attention to the physical part of the correlator functions (6) and (7). Writing 
a dispersion relation in the variable (p + q) 2 , one can separate the B meson ground state 
contribution to the correlator functions il^ 1 -* and H^ 2 \ by inserting a complete set of states 
between the currents in (6) and (7) focusing on the term \B) (B\ : 
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(22) 
(23) 



Using the definitions of the formfactors (see eqs.(3) and (4)) in (22) and (23) and equating 
these expressions to eqs.(19) and (20), we get the sum rules for the formfactors. The 
remaining part of the calculation follows from the QCD sum rules procedure: perform the 
Borel transformation on the variable (p + q) 2 and subtract the continuum and higher states 
contributions invoking quark-hadron duality. (Details of these procedures can be found in 
[17-19]). After this procedure we obtain the following sum rules for the formfactors: 
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From eq.(26) and eq.(27) we get a new relation between formfactors A 3 , A and A 2 : 

a ( 2n a t i\ A 2 (q 2 )q 2 

A 3 {q - A (q = - ■ . 

2m k* [tub + mK*) 

For the formfactors 7\, T%, and T 3 , we get the following sum rules: 
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Using the equation of motion we can relate T3 and A3 — Aq by: 



T 3 (q 2 ) = m K *(m b 



m s 



A 3 (q 2 )-A (q 2 



(32) 



Here M is the Borel mass parameter. The lower integration limit 5 



depends on the 



effective threshold So above which the contributions from higher states to the dispersion 
relation (22) and (23) are cancelled against the corresponding piece in the QCD represen- 
tation (19) and (20) Note that the sum rules for V(q 2 ) and Ai(q 2 ) and Ti(q 2 ) in the light 

Our results agree with that of fl7fl. The region of appli- 



cone QCD are derived in [17 



cability of these sum rules is restricted by the requirement that the value of q 



mi 



be 



sufficiently less than zero. In order not to introduce an additional scale, we require that 
q 2 — ml < (p + q) 2 — m b which translates to the condition that m\ — q 2 is of the order of the 
typical Borel parameter M 2 ~ 5 4- 8 GeV 2 . From this condition we obtain that the region 
of applicability of the sum rules is q 2 < 15 -4 17 GeV 2 , which is few GeV 2 below the zero 
recoil point. 
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Finally we calculate the differential decay rate with longitudinal polarization of the final 
leptons. The differential decay rate is given by: 
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where A = 1 + r 2 + s 2 — 2r — 2s — 2rs, r = s = -4-, £ is the longitudinal polarization 



of the final lepton, mi and v = \J 1 — are its mass and velocity, respectively. In eq. (33) 
A, Bi, B 2 , B 3 , C, Di, D 2 , and D 3 are defined as follows: 
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D 3 = -C W ^(A 3 - Ao) . 

For the dileptonic decays of the B mesons, the longitudinal polarization asymmetry, Pl, of 
the final lepton, /, is defined as 



P L{f) = g K """^ 11 , (34) 

where £ = — 1(+1) corresponds to the left (right) handed lepton in the final state. In the 
Standard Model, this polarization asymmetry comes from the interference of the vector or 
magnetic moment and axial vector operators. If in eq.(33) the lepton mass is equated to 
zero, our results coincide with the results in and if mi ^ they coincide with the results 
in IBH. 



3 Numerical Analysis 

For the input parameters which enter the sum rules for the formfactors and the expressions 
of the decay width we have used the following values : 

m b = 4.8 GeV, m c = 1.35 GeV, m„ = 0.106 GeV, m T = 1.78 GeV, 



A Q cd = 225 MeV, m B = 5.28 GeV, m K » = 0.892 GeV, s = 36 GeV 2 , M 2 = 8 GeV 2 

In Fig.l we present the q 2 dependence of the formfactors V(q 2 ), Ai(q 2 ), A 2 (q 2 ) and 
A (q 2 ) (the formfactor A 3 can be easily obtained from eq.(28)). All these formfactors 
increase with q 2 . From these figures we see that A 2 (q 2 ) increases with q 2 strongly, but 
Ai(q 2 ) and A (q 2 ) do so smoothly. At this point let us compare our results on these 
formfactors with the results which are obtained from 3-point QCD sum rules analysis in 
p. In our case Ai(q 2 ) increases with q 2 , but in |J it decreases with q 2 . The behaviour of 
the other formfactors are similar. 

In Fig. 2 we depict the dependence of the formfactors 7\, T 2 , and T 3 on q 2 . In this 
case also all formfactors increase with q 2 . For formfactors T 2 and T 3 , our predictions 
on their q 2 dependence also differ from the predictions of 0. In ||, T 2 is positive and 
smoothly decreases, the value of T3 is negative for all q 2 . Note that our predictions on the 
q 2 dependence of all formfactors coincide with relativistic quark model predictions |TTJ . The 
source of discrepeancy of our results with the predictions of on Ai, T 2 , and T 3 should 
be carefully analysed. This lies out of the scope of this paper. We are planning to come 
back to the analysis of these points in our forthcoming works. 

In Fig.3(4) we present the q 2 dependence of the branching ratios for B — > K*fi + fi~ 
(B — > K*t + t~) decay, with and without the long distance effects, respectively. 

In Fig. 5 we plot the longitudinal polarization asymmetry Pi as a function of q 2 for 
B — > K* and B — > K*t + t~ , with m t = 176 GeV, with and without the long distance 



9 



effects. From this figure we see that Pi vanishes at the threshold due to the kinematical 
factor v and that the value of Pi for B — > K* ji + pT decay varies in the region (—0.5 + +0.5), 
when the resonance ip, if)' mass region is excluded. In the B — ► K*t + t~ decay case, without 
long distance effects Pi is negative for all values of q 2 , and only in the resonance ip' mass 
region Pi become positive. Therefore the study of the longitudinal polarization Pi can be 
very useful for understanding the relative roles of the long and short distance contributions 
in the B — > K*l + l~ decay. 

At the end of this section let us compare our results on the Branching ratio of the 
decay with those in [§, 11 1. The value of the branching ratio is close to the 

In our opinion this is due to the 



T 



B -> K*l 
results of ^ 
over estimation of the formfactors in Wl 



but about 30 times smaller than that of ITT 



4 Conclusions 

We calculate the transition formfactors for the exclusive B — > K*l + l~(l = fi, r) decay in the 
framework of the lightcone QCD sum rules, and investigate the longitudinal polarization 
asymmetries of the muon and tau in this decay. It is shown that some of the formfactors 
in light cone and 3-point QCD sum rules have absolutely different q 2 dependence. It is 
found that the value of the longitudinal polarization Pi(fi) in the region (—0.5; +0.5) and 
Pl(t) in (0; —0.6). We also calculate the integral branching ratios and find that they are 
Br(B -»■ K*fi + fi~) = 12.06 and Br(B -> K*t + t-) = 0.217. 

Few words about the possibility of the experimental observation of this decay are in 
order. Experimentally, to observe an asymmetry Pi of a decay with the branching ratio 
Br at the na level, the required number of events is iV = B ™ pi (see [0). For example, to 

observe the r lepton polarization at the exclusive channel B — > K*t + t~ at the 3a level, one 
needs at least N = 1.66 x 10 9 B B decays. Since in the future S-factories, it is expected 
that ~ 10 9 B-mesons would be created per year, it is possible to measure the longitudinal 
polarization asymmetry of the r lepton. 
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Figure 3 
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Figure 5 
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Figure Captions 

1. The q 2 dependence of the formfactors V(q 2 ), Ai(q 2 ), A 2 {q 2 ) and A (q 2 ). 

2. The q 2 dependence of the formfactors Ti(q 2 ), T 2 (q 2 ) and T 3 (g 2 ). 

3. a) Invariant mass squared distribution of the lepton pair for the decay B — > K*fi + fi~ 
which includes only the short distance contributions. 

b) The same as in a) but including long distance effects, too. 

4. The same as in Fig. 3, but for B — > K*t + t~ decay. 

5. a) The longitudinal polarization asymmetry P L for the B — > K*fi + fi~ decay, 
b) The same as in a), but for B — > K*t + t~ decay. 

In these figures, solid line corresponds to the short distance contributions only and dashed 
to the sum of both short and long distance contributions. 
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